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The maximum solubility of titanium in BeO is estimated on the basis of calculations of the fonnation enthalpy 
of beryllium oxide, doped with titanium, performed using quantum-chemical calculations within the frame¬ 
work of the DFT theory. A method for estimating numerically the maximum solubility of impurities in crystals 
is examined for beryllium oxide doped with titanium. The method is based on using the results of non- 
empirical quantum-chemical calculations. The method can be used to estimate the maximum solubility of di¬ 
verse impurities in inorganic compounds and ceramic materials based on these compounds. 


Ceramic based on beryllium oxide BeO possesses unique 
physical-chemical properties, combining high chemical, 
thermal, and radiation resistance, thermal conductivity, and a 
number of other characteristics. This makes it a promising 
material for use in electronics, new areas of engineering, and 
special instrument building [1]. 

It is known [1] that beryllium oxide essentially does not 
form solid substitution solutions. According to published 
data [1,2], the solubility of different elements fluctuates in 
the range 10" 3 — 10^ 4 % (atomic content). On the other 
hand, impurity defects can have a large effect on the phy¬ 
sical-chemical properties and therefore also on the functional 
characteristics of beryllium oxide (luminescence, exoemis¬ 
sion, electric conductivity, and others) [1]. Consequently, it is 
of interest to develop methods and approaches which make it 
possible to estimate the maximum solubility of various ele¬ 
ments in beryllium oxide. 

We have examined the possibility of using the results ob¬ 
tained within the framework of modem quantum-chemical 
computational schemes based on the electronic density func¬ 
tional theory (DFT) to estimate the maximum solubility of 
impurities in BeO for the example of the system beryllium 
oxide - titanium impurity (BeO - Ti). 

The formation in beryllium oxide of an impurity defect 
caused by inclusion of titanium Ti Be was modeled by substi¬ 
tuting a titanium atom for a beryllium atom in the initial 
crystal. In turn, the beryllium atom is removed from the be¬ 
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ryllium oxide lattice. This can be represented by the formal 
reaction 

BeO + xTi -A Bej _ x Ti v O + xBe. (1) 

For subsequent energy estimates, we performed a quan¬ 
tum-chemical calculation of the electronic structure of the re¬ 
gions indicated above: BeO and Be, _ x Ti x O. We used 
so-called supercells, which model defect-free beryllium 
oxide and BeO with a point defect Ti Be (Be] x Ti v O). Their 
energy spectrum was calculated ah initio by the VASP-PAW 
method [3] with a nonlocal generalized gradient approxima¬ 
tion of the exchange-correlation potential (GGA) in the PBE 
form [4]. The method of tetrahedra was used to integrate 
over the Brillouin zone [5]. The energy of the plane waves 
was 400 eV, which gave good convergence. The relaxation 
effects which appear when an impurity atom is embedded 
were taken into account by calculating the forces F acting on 
the atoms. The convergence criterion when optimizing the 
geometric parameters of a cell corresponded to the condition 
F < 0.005 eV/A. 

To determine the optimal size of the supercells, calcula¬ 
tions were performed for a series of BeO cells with wurtzite 
structure (B4, spac group C6mc), containing 16, 32, 72, and 
128 atoms. The structure B4-BeO was given by three pa¬ 
rameters: the lattice constants a and c and an internal para¬ 
meter u which determines the relative displacement of the 
sublattices of the O and Be atoms along the c axis. The total 
energy E tol of these supercells as a function of the number of 
atoms is displayed in Fig. 1. It was found that E, ot changes 
very little as the number of atoms increases above 72. Conse- 
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Fig. 1. Total energy E tot of a super cell of wurzite-like beryllium 
oxide versus the number of atoms. 


quently, 72-atom BeO - Be 36 0 36 supercells were used in 
subsequent calculations. The solid solution Be, _ z Ti x O, sim¬ 
ulated by replacing a beryllium atom with a titanium atom, 
for such a cell (Be 35 Ti 36 ) corresponded to the formal 
stoichiometry Be 0 972 Ti 0 028 O. 

As a result, the values of the total energies of 
Be 0 972 T 10 02 sO and BeO were obtained and found to be 
-511.639720 and- 512.615290 eV, respectively. 

The formation enthalpy h Ae{ of a point defect Ti Be in the 
system BeO - Ti was estimated as 

Kef = E m (Be 35 Ti0 36 ) - E tol (BeO) 36 + A//(Be) - AH (Ti), 

where E tot (Be 35 Ti0 36 ) and E tot (BeO) 36 are the total energies 
of the supercells simulating Be 0 972 Ti 0 028 O and BeO, respec¬ 
tively; AH (Be) and AH (Ti) are the heats of formation of be¬ 
ryllium and titanium at the calcination temperatures of the 
ceramic. 

These quantities were taken from the reference literature 
[5], Then h def ® 84 kJ/mole for T « 1500 K. 

The Gibbs formation energy of a single defect g def can be 
determined from the relation 

*^def ~ ^def - ^^def ’ 

where AS def is the entropy change due to the formation of a 
defect. 

In turn, the quantity A5 def can be estimated from the 
equation 

Av 0 
i _ fcT 

AS de{ = A S l + AA 2 = 3 R In ~ e — + [S (Be) - 5(Ti)], (2) 
1-e 

where R is the universal gas constant, h is Plank’s constant, 
and k is Boltzmann’s constant. 

The first term AS) takes account of the entropy change in 
Einstein’s approximation, due to the difference of the fre¬ 
quencies v 0 and v, characterizing the vibrational motion in a 
crystal without and with defects, respectively. According to 


the data in [2], v 0 « 1500 cm 1 and v ~ 1000 cm 1 . The cal¬ 
culations showed that AS) « 5 J/(mole • K). 

The second term A S 2 in equation (2) takes account of the 
difference of the entropies of the final (Be) and initial (Ti) 
states according to the reaction (1). This quantity can be esti¬ 
mated from the thermodynamic reference data [5], according 
to which AS 2 * - 35 J/(mole ■ K) at T~ 1500 K. 

Thus, the estimates made for a point defect Ti Be in the 
system BeO - Ti lead to the values AS def « - 30 J/(mole • K) 
andgdef * 130 kJ/mole. 

When a number of defects Ti Be form, the change AG in 
the Gibbs energy of the system BeO - Ti will be given by the 
relation 


^defSdef ^>conf’ 


where N d f is the number of defects in the crystal and AS def is 
the change in the configuration entropy of the crystal due to 
the formation of defects. 

The quantity AS conf is determined by the relation [6, 7] 


ASconf^l* 


N ! 


N ief \(N-N det )\ 


( 3 ) 


where N is the total number of sites in the crystal. 

The equilibrium value /V^ ( . is obtained from the condi¬ 
tion 


d(AG) 

dN det 


( 4 ) 


The solution of equation (4) obtained using Stirling’s for¬ 
mula in the relation (3) has the form 


_ * £ 
1 N 


= exp 


& def 

RT 


where c at is the equilibrium atomic fraction of impurity de¬ 
fects Ti Be in a BeO crystal. 

As a result, for the defects Ti Be at T~ 1500 K, taking ac¬ 
count of the estimates presented, g'def c at S! 3 • 10 ~ 5 , which 
agrees well with published data on the low solubility of im¬ 
purities in the beryllium oxide lattice [1,2], 

In summary, a method for obtaining a numerical estimate 
of the maximum solubility of impurities in crystals, which is 
based on using the results of nonempirical quantum-chemi¬ 
cal calculations, was examined for beryllium oxide doped 
with titanium. The method can be used to estimate the maxi¬ 
mum solubility of diverse impurities (for example, d and / 
elements) and inorganic compounds and ceramic materials 
based on them. 
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